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• . We analyze the correlation of the positions of 7-ray sources in the Fermi Large 

Q^. Area Telescope First Source Catalog (IFGL) and the First LAT Active Galactic 

Q . Nuclei (AGN) Catalog (ILAC) with the arrival directions of ultra-high-energy 

cosmic rays (UHECRs) observed with the Pierre Auger Observatory, in order to 
c^ . investigate the origin of UHECRs. We find that Galactic sources and blazars 

identified in the IFGL are not significantly correlated with UHECRs, while the 
ILAC sources display a mild correlation (2.6cr level) on a ~ 2.4° angular scale. 
t~^ . When selecting only the ILAC AGNs closer than 200 Mpc, we find a strong 

ff-^ . association (5.4(t) between their positions and the directions of UHECRs on a 

Ri 17° angular scale; the probability of the observed configuration being due to 
an isotropic flux of cosmic rays is 5 x 10~^. There is also a 5cr correlation with 
O '. nearby ILAC sources on a 6.5° scale. We identify 7 "7-ray loud" AGNs which are 

associated with UHECRs within ^ 17° and are likely candidates for the produc- 
tion sites of UHECRs: Centaurus A, NGC 4945, ESQ 323-G77, 4C+04.77, NGC 
1218, RX J0008.0+1450 and NGC 253. We interpret these results as providing 
additional support to the hypothesis of the origin of UHECRs in nearby extra- 
galactic objects. As the angular scales of the correlations are large, we discuss 
the possibility that intervening magnetic fields might be considerably deflecting 
the trajectories of the particles on their way to Earth. 

Subject headings: gamma rays: observations — cosmic rays — galaxies: active 
— galaxies: jets — methods: statistical 
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Introduction 



The nature of the ultra-high-energy cosmic rays (UHECRs) with energies abov e ^10 



20 



eV re mains enigmatic since they were first observed, more than half a century ago (JLinsley 
19631 ). If such cosmic rays are com posed predom i nantly of protons or nucle i, the Greisen, 
Zatsepin and Kuzmin (GZK) effect (JGreisenlll966l : IZatsepin &: KuzminI 119661 ) restricts their 
possible sources to the nearby unive rse, closer than about ~ 100 — 200 Mpc from Earth (the 
"GZK horizon", iHarari et al.ll2006l ). Under the assumption that the UHECRs sources are 
relatively nearby and not uniformly distributed on the sky, then an anisotropic distribution 
of arrival directions is expected, as long as the deflections caused by the intervening magnetic 
fields are sufficiently small. 



W ith the unprecedented capabilities of the Pierre Auger Observatory (PAO: I Abraham et al, 
2004r). it is now possi ble to map the arrival directions of UHECRs with a ~ 1° precision 
( lAbraham et al.l 120071 ). This makes the cross-correlation analysis of the arrival directions 



with catalogs of astronomi c al objects a powerful tool for tracking the UHECR production 
sites (lAbraham et al.l 120071 : lAublinl l2009l : lHaguell2009l ). For instance, the anisotropy in the 
arrival directions of UHECRs was recently demonstrated with a significance level above 99%, 
through the correlation between the directions of such particles and t he positions of nearby 
activ e galactic nuclei (AGNs) from the Veron-Cetty and Veron catalog (jAbraham et al.ll2007l . 
20081 ). This indicates that the highest-energy particles have an extragalactic origin within the 
GZK horizon, with either AGNs or other objects with a similar spatial distribution being the 
likely astrophysical sources of UHECRs. Note that the High Resolution Fly's Eye (HiRes) 



exper iment does not confirm the correlations suggested by the PAO data (jAbbasi et al.ll2008 
2010h . 



Since the pioneering work of lAbraham et al.l ( 120071 ). several authors have been studying 
the correlation of the PAO events' directions with different classes of astrophysical objects. 
There is circumstantial evidence tha t at least some UHECR s appear to be associate d with 
relatively few nearby radio galaxies ( jNagar fc Matulichll2008t iMoskalenko et al.ll2009l). with 
a puz zling paucity of UHECRs in the direction of the Virgo cluster ( IZaw. Farrar fc Greene 
20091 ). There is also evidence that spiral g alaxies seem to host the producers of ultra-high 



energy cosmic rays (JGuisellini et al 



20081 ) . IMoskalenko et al.l ( 120091 ) noted the importance 



of taking into account the AGN morphology in correlation studies. Moskalenko et al. also 
pointed out the possibility of UHECRs being correlated with la rger defiection angles and/or 
more distant sources than considered by lAbraham et al.l ( 120071 ). An understanding of which 
particular class of extragalactic sources is associated with the production sites of UHECRs 
is clearly missing. 



Astrophysical 7-ray sources are appealing candidate accelerators of UHECRs, since by 
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their very nature they are non thermal objects where lar ge energy transfers and extreme 
particle acceleration take place (JDermer &: Razzaquell2010l ). Thanks to the Large Area Tele- 
scope (LAT) aboard the Fermi Gamma-ray Space Telescope [Fermi; lAtwood et al.l |2009[ ) 
which has been conducting a sky survey at ene r gies ab ove 100 MeV, we have an unprece- 
dented map of the 7-ray sky (JAbdo et al.ll2009al . l2010al ). In its first 11 months of operation 
Fermi detected and characteri zed 1451 7-ray obj ects in the sky, culminating in the Fermi 
LAT First Source List flFGLi lAbdo et al.l[2010ah and the First LAT AGN Catalog (ILAC; 
Abdo et al.ll2010bl ). These detailed maps of the 7-ray sky are a rich dataset which that be 
used in cross- correlation analyses to shed light on the origin of energetic cosmic rays. 



Mirabal &: Oyal (120101 ) recently invest igated the correlation of the complete sample of 
IFGL sources with the PAO dataset (as of lAbraham et al.ll2008l ). Their analysis was carried 
out without any redshift cut-off or type selection on the IFGL, and counted as sociations when 
PAO events occur within circles of radius 3.1° around the Fermi sources. iMirabal fc Oya 
(I2OIOI ) found no evidence that UHECRs are associated with Fermi sources for the particular 
angular radius considered. This does not imply though that there is no association at all 
between 7-ray sources and UHECRs. Considering that the IFGL catalog is a "mixed bag" 
of 7-ray sources with quite different properties, if one particular class of 7-ray emitters is 
responsible for producing UHECRs, then other sources in the catalog bearing no relation 
with the cosmic rays might introduce "noise" in the correlation analysis. Furthermore, the 
correlation can be weak at a given angular separation and stronger at others. 

In this work we use the IFGL and ILAC catalogs produced by the Fermi LAT to 
investigate the correlation between the positions of 7-ray sources and the arrival directions 
of UHECRs measured by the PAO. Our goal is to unveil if specific types of 7-ray sources 
might be driving any correlation and at what angular separations these correlations become 
significant. In ^we describe the datasets we use. In ^we describe our cross-correlation 
method. In ^we list the results of our cross-correlation analysis when applied to the different 
subsamples of the IFGL and ILAC catalogs. We discuss the implications of our results in 
^ In particular, we list the potential "7-ray loud" UHECR accelerators and discuss their 
properties in §5.11 Section [6] closes with our concluding remarks. 



Data 



We base our analysis on the distribution of arrival directions of UHECRs with energies 
exceeding 5.7 x 10^^ eV, collected by the surface array of PAO between 1 January, 2004 
and 31 Aug ust 2007, with an integrated exposure of 9.0 x 10^ km sr year ( Abraham et al. 
20071 . I2OO8I ). This data set corresponds to 27 events with zenith angles smaller than 60° and 
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angular resolution of ~ 1° fjAvd 120071 ) . 



To search for the pot ential astrophysical sources of UHECRs, we use the Fermi LAT 
First Source List (IFGL; kbdo et al.l boiOaf ) and the First LAT AGN Catalog (ILAC; 
Abdo et al.l l2010bl ) of 7-ray sources produced after the first eleven months of operation of 
Fermi. Source detection is based on the average flux over the ll-month period corresponding 
to a statistical significance higher than 4cr. The IFGL consists of the 1451 sources detected 
and characterized in the 100 MeV to 100 GeV range. A subset of the IFGL consisting of 
671 sources constitutes the ILAC. 



1037 sources in the IFGL are within the PAO field of view. lAbdo et al.l (l2010al ) pro- 
vides identifications or plausible associations of these 7-ray sources with objects in other 
astronomical catalogs: 93 of them correspond to Galactic sources (pulsars, pulsar wind neb- 
ulae, supernova remnants. X-ray binaries and globular clusters) and 453 to AGNs in the 
ILAC. Of the ILAC AGNs, 367 are blazars, 69 are AGNs of uncertain type and 17 are non- 
blazar AGNs that i nclude, for example, the radio galaxies M87, Centaurus A and 3C 207.0 
(JAbdo et al.ll2010bl ). 490 IFGL sources could not be associated with any counterpart. Most 
of the unassociated sources are near the Galactic plane, and c ould not be identified beca use 
of a combination of Galactic extinction and source confusion ( Abdo et al.ll2009al . l2010a| ]bl). 



Figure [T] shows the sky map in Galactic coordinates of the UHECRs events together 
with the IFGL a nd ILAC sources. Also shown is the supergalactic plane along which nearby 
galaxies cluster (JLahav et al.ll2000l ). 



Cross-correlation method 



In order to quantify the statistical cross-correlation between the Fermi and PAO data, 
we compute the cumulative number of associated cosmic ray events as a function of angular 
separation tp (hereafter CCF for cumulative correlation function) from the data sets. The 
CCFs are computed by drawing circles of increasing radius around the Fermi sources, and 
counting the number of cosr nic ray events t h at occ ur within such circles. This is qualitatively 
similar to the approach of [Stephen et al.l (120051 ) for the correlation analysis between two 
different X-ray catalogs. 

The significance of the measured CCF (number of correlations measured in the observed 
data No as a function ip) is evaluated by comparing it against that expected by chance cor- 
relations assuming an isotropic flux of UHECRs. The chance CCF - number of correlations 
expected by chance A''^ vs. ^ - is obtained by building simulated sets of UHECRs with the 
same number of observed events, drawn from an isotropic distribution of arrival directions 
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convolved with the exposure function of the PAO (jSoniniersll200ll ) 



For a given configuration of Fermi sources and a given ■?/', our null hypothesis is that the 
number of associations measured in the data is the same as the average number expected by 
chance: No = (Nc). We calculate the p- value P as the probability that Nc exceeds No- In 
practice, this p-value can be interpreted as the likelihood that the associations measured in 
the data could have occurred by chance if the null hypothesis is true, i.e. the likelihood that 
the observed configuration of 7-ray sources and UHECRs is a random accident. Assuming 
that the distribution of N^ values obtained from the Monte Carlo simulation is normal, we 

define P as 

2 [1 - <l>(AA^/a)] , AN>0 



^^' { 2$(AiV/a), AN <0 

where AN{tp) = N^, — (Nc), cr(^) is the standard deviation of the distribution of values of 
Nc and $(a;) is the cumulative distribution function for a normal distribution. 

When describing our results in the sections that follow, we will often refer to the sta- 
tistical significance level with which we can reject the null hypothesis. We will refer to this 
significance level as S* = v2 erf^ (1 — P), given in standard deviations. S is just another 
way of expressing the p-value P. 



3.1. Consistency check 

As a consistency check of the method outlined above, we applied it to calculate the 
cross-correlation between the arrival directions of UHECRs and the 292 AGNs with redshift 
z < 0.017 within the PAO field of view listed in t he 12th edition of the Veron-Cetty and 



Veron (VCy) catalog (IVeron-Cetty fc Veronll2006l ). previously studied by [Abraham et al. 



(I2OO7I . I2OO8I ). Figure [2] shows the resulting correlation signal. The upper panel in Figure [2] 
shows the measured CCFs, with angular separations ranging from 1° (approximately the la 
uncertainty in the arrival directions of PAO cosmic rays) up to 10°. The lower panels display 
the probability P. 

Figure [2] illustrates that P reaches a minimal value of 8.3 x 10~^^ at the angular sepa- 
ration ^ 3.2°, with 20 cosmic ray events correlated with AGNs at this angular scale, while 
only 5.9 correlations are expected by chance. Therefore, the results of our consistency check 
confirm with at least a significance level S = 6.5a the anisotropy of the arrival directions of 
UHECRs and are in excellent agree ment with those obtained by the statistical method used 



by the Pierre Auger Collaboration (lAbraham et al.l 120071 . 120081 ). This lends strong support 



;ai nea py 
)7l. l2008h . 



to the validity of our statistical method to probe the degree of correlation between the PAO 
and Fermi data sets. 
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As a matter of fact, when we apply our statistical test to the revised AGN sample in the 
13th edition of the VCV catalog using the same redshift cutoff and considering only AGNs 
within the PAO field of view, we find that the association with UHECRs is maximized at 
tl) K. 2.3° with Pmin = 6.7x 10"^*^ and correspondingly S'max = 6.2(j. Therefore, the correlation 
signal slightly changes when cross-correlating PAO events with the updated VCV catalog. 



4. Results 

We applied the above method to obtain the cross-correlation signal between the IFGL/ILAC 
and PAO sources. In the following subsections we analyse the correlation signal of PAO 
events with: (1) all sources in the IFGL catalog ( §4.ip . (2) only to IFGL Galactic sources 
(S21) and (3) only to ILAC AGNs (Sections |0] and IMD- In S3]we impose no restriction 
on the redshifts of the sources while in §4.4l we restrict ourselves only to nearby AGNs within 
the GZK horizon. 



4.1. UHECRs vs. all IFGL sources 

Figure [3] shows the correlation signal resulting from comparing the positions of all 1037 
IFGL 7-ray sources within the PAO field of view with the directions of energetic cosmic ray 
events. The upper panel shows the measured CCFs, with angular separations ranging from 
1° (approximately the la uncertainty in the arrival directions of PAO cosmic rays) up to 
10°. The lower panels display the probability P. 

Both panels of Figure |3] show that the number of observed correlated events follows the 
pattern of isotropic expectations within the la significance until ip k. 4° . In the range of 
angular separations 4° < ^ < 6°, there is a systematic paucity of measured correlations with 
respect to chance expectations with 5* > la. In particular, the maximal "anti-association" 
between the IFGL and PAO datasets corresponds to the global minimum of the probability 
signal Pmin = 0.6%, which occurs aX ip ^ 5.1°. This minimal probability corresponds to the 
significance level S = 2.7a. At ip ^ 5°, 22 PAO events correlate with all IFGL sources while 
25.8 are expected from an isotropic flux. 

The case corresponding to circles of radius 3.1° deserves special attention. This is 
the estimated angular radius around AGNs in the Veron-Cetty ar id Veron catalog tha t 



maximizes their correlation with the arrival directions of UHECRs (jAbraham et al.l 120071 ). 



Figure E] illustrates that at ip = 3.1° we find a high probability of ch ance correlation, P = 0.6 



{S = 0.5a). This agrees with the analysis of iMirabal fc Oyal (120101 ). 
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The IFGL catalog consists of a "mixed bag" of Galactic and extragalactic 7-ray sources 
with quite different properties. In the subsections below, we iii v estigat e the correlation of 



the individual classes of 7-ray sources identified by lAbdo et al.l (l2010al ) with the positions 
of PAO events. 



4.2. UHECRs vs. IFGL Galactic sources 

FigureHlshows the correlation signal measured in the subset of 93 IFGL Galactic sources. 
The number of observed correlated events follows well the pattern of an isotropic expectation 
within the la significance. There is a narrow fluctuation in the signal above S* = la at 
if} K. 2.4°. Nevertheless, there are only 3 correlated events at this angular distance with 1 
correlation expected from isotropy. 



4.3. UHECRs vs. ILAC AGNs 

In this section we study the likelihood of association between UHECRs and the 453 
AGNs identified in the IFGL, i.e. the sources in the ILAC catalog within the PAO field of 
view. Figure |5] shows the sky map in Galactic coordinates of the PAO events and the ILAC 
sources. These AGNs are located at distances ranging from z ?» 0.001 (the nearest ILAC 
AGN is NGC 253) vcp to z ^ 3.2, which corresponds to the FSRQ blazar PKS 0336-017 



(lAbdo et al.ll2010bl ). 



The left panel of Figure |6] shows the correlation signal of the PAO events with all ILAC 
AGNs. There is a maximal association between the two data sets at the angular scale 
ip ^ 2.4° (i.e. P reaches a minimal value Pmin ~ 0.01, corresponding to ^max = 2.6a). 
Within 2.4°, 12 events correlate with all the Fermi AGNs without distinction of distance or 
AGN class, whereas 6.5 correlations are expected by chance. 

The middle panel of Figure |6] shows the correlation signal with the 367 ILAC blazars. 
There is no significant correlation (i.e. S < la) between blazars and UHECRs that cannot 
be explained by an isotropic flux of cosmic rays. We also study the amount of correlation 
obtained when cross-correlating separately the two classes of blazars, BL Lacertae and FSRQ, 
with the PAO UHECRs. We obtained for the 169 BL Lacertae ILAC sources a quite similar 
correlation signal to that of all blazars (i.e. consistent with chance correlations). If we select 
only FSRQ sources, we have a sample of 198 FSRQ objects, which increases the significance 
of the correlation up to 2.3a at t/) ~ 6.5°. It is interesting to note that ILAC FSRQ blazars 
are located at ^ > 0.15, well outside the GZK horizon. Therefore, we should not expect an 



association of these objects with the arrival directions of UHECRs. 

The right panel of Figure E] displays the amount of cross- correlation between the subset 
of 86 AGNs consisting of non blazars (i.e. 17 misaligned jet sources such as the radio 
galaxies Centaurus A and M87) together with the AGNs that have not been identified yet 
(69 sources). The lower panel shows that the correlation significance raises above la in 
the range 2° < t/^ < 7.5°. In particular, the probability signal reaches the minimum value 
-Pmin ~ 5 X 10""^ (5'max = 3.5(t) at z/' ~ 3.2°, with 8 observed PAO event correlations, whereas 
only 2.9 chance correlations are expected. These PAO events are correlated with 7 AGNs, of 
which Gen A, 4G+04.77 and NGG 4945 are the only sources whose AGN types are known. 



4.4. UHECRs vs. nearby ILAC AGNs 

Protons or nuclei with energies above 60 EeV interact with the cosmic microwave back- 
ground and suffer the GZK effect : a st rong attenuation of their flux for distant sources 
(lGreiserull966l : IZatsepin fc Kuzmirull966f) . This defines a GZK horizon located roughly at a 
distance ~ 100 — 200 Mpc (JHarari et al.ll2006l ). Astrophysical sources located outside the 
GZK horizon are not expected to contribute appreciably to the observed flux of energetic 
cosmic rays. Motivated by this expectation, we investigate the correlation of PAO events 
with ILAC AGNs located within the GZK horizon. We define the GZK horizon as being 
located at a comoving distance of 200 Mpc, which corresponds to z^^^^ = 0.048, assuming 
the AGDM cosmology with Hq = 71 km s'^ Mpc"\ fi^ = 0.270, ^a = 0.730 and fi^ = 
jKomatsu et al.ll2009h . 



We use the redshift information provided in the ILAC (lAbdo et al.ll2010bl ) to discard 
the AGNs within the PAO field of view for which z > Zmax- Only 9 AGNs remain, including 
the radio galaxies Centaurus A (the nearest radio galaxy) and M87, the Seyfert 2 NGC 
49 45, the Seyfert 1 ESQ 323- G77 and the starburst galaxy NGC 253 (as previously noted 
by lDermer &: Razzaqudl2010l ). Only three AGNs within the GZK horizon are BL Lacertae 
objects. There is no FSRQ AGN within the GZK horizon, with the nearest FSRQ blazar in 
the ILAC being located aX z = 0.15. ILAC sources inside the GZK horizon are shown as 
squares in Figure [51 

Figure [7] shows the correlation signal between ILAC AGNs with z < 0.048 and UHECRs, 
with ip reaching up to 20°. The noteworthy feature here is that there are three characteristic 
angular scales where the probability of the observed configuration being explained by an 
isotropic flux of cosmic rays is minimized: ipi ~ 2.3°, ^2 ~ 6.5° and ip^ ~ 16.9°. The 
first local minimum of the probability signal located at ipi ^ 2.3° corresponds to only 4 
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associations, with P{ipi) = 6 x 10^^ and S{il)i) = 4a. The second local minimum of the 
probability signal located at iIj2 ~ 6.5° corresponds to 7 measured correlations whereas 0.9 
chance correlations are expected, with P{ip2) = 5 x 10~^ and 5(^2) = 5cr. The absolute 
minimum of P is characterized by ip^ ^ 16.9°, No = 17, N^ = 5.7, P{ipz) = 5 x 10^^ and 
S{^^) = 5.4a. 

Within the distance ipi, 4 PAO events correlate with 3 AGNs: NGC 4945, Centaurus 
A and 4C+04.77; at ^^2, 7 PAO events correlate with 4 AGNs: NGC 4945, Centaurus A, 
ESQ 323-G77 and 4C+04.77; at V's, 17 events are associated with 7 AGNs: NGC 4945, 
Centaurus A, ESQ 323-G77, 4C+04.77, NGC 253, NGC 1218 and RXJ0008.0+1450. Table 
[1] lists the 9 ILAC AGNs located inside the GZK horizon, including their IFGL names. 
Galactic co o rdinate s (longitude / and latitude 6), identifications and redshifts according to 
Abdo et al.l ( 2010al Jbl). We also list the number of PAO events potentially correlated with 



each AGN at the angular scales ip2 and ip^. The ellipses in Figure |5] correspond to circles of 
radius ^/'s centered on the GZK AGNs. 



Discussion 



We begin by pointing out some caveats in our analysis of the correlation of ILAC 
AGNs and PAO events. Fi rstly, the ILAC ca talog of AGNs is increasingly incomplete 
towards the Galactic plane (jAbdo et al.l l2010al Jbl: see Figure [5]) due to a combination of 
the effects of dust extinction in the Milky Way and source confusion. These effects may 
have some impact on our estimate of the correlation strength, since we might be missing 
some AGNs lying close to the Galactic plane. Secondly, there are many AGNs in the ILAC 
that have not ye t been identified w ith counterparts in other wavelengths, and lack redshift 
determinations (jAbdo et al.l l2010bl ). This may affect our cross-correlation analysis using 
the subset of ILAC GZK AGNs, since we might be missing AGNs inside the GZK horizon 
that are yet unidentified in the ILAC. With the progressive increase in the exposure time 
of the Fermi LAT survey observations, future releases of the Fermi LAT source catalog 
will improve on both these issues: they will increase the number of AGN detections near the 
Galactic plan e and increase the number of associations of 7-ray sources with multiwavelength 
counterparts JAbdo et al.ll2009a[ l201o"al Jb[). 



It is interesting to point out that ip = 'i.2° - the angular separation which maximizes the 
association of PAO events with the AGNs in the 12th VCV catalog - does not correspond 
to a local minimum of the probability signal in our analysis of the correlation of ILAC GZK 
sources and PAO events, even though at 3.2° the hypothesis of the observed correlations 
being due to chance can be rejected at the 4.2cr level. On the other hand, the radius 2.3° 
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that maximizes the association with the sources in the 13th VCV catalog corresponds to a 
local minimum of the correlation signal in Figure [3 

We find no correlation between ILAC blazars and UHECRs when considering all blazars 
detected by Fermi/LAT up to ^ ~ 3.2. In other words, relativistic jets nearly pointed along 
our line of sight do not seem to be significantly associated with the arrival directions of 
UHECRs on any angular scale. This result is in disagreement with the controversial claims 
of correlation of HiR es data with BL Lac object s (JGorbunov et al.ll2004J ). but is expected on 
theoretical grounds (Ivan Putten fc Guptal 120091 ). 



5.1. The potential UHECR accelerators 

Our results pinpoint 7 "7-ray loud" ILAC AGNs as the potential sites for the production 
of PAO UHECRs (see Tabled]): Cen A, NGC 4945, ESQ 323-G77, 4C+04.77, NGC 1218, 
RX J0008. 0+1450 and NGC 253 (ordered by increasing number of associated UHECRs at 
the angul ar distance 16.9°). N one of these AGNs display significant variability in their 7-ray 
emission kbdo et al.l boiOal Jbll The three sources Cen A, NGC 4945 and ESQ 323-G77 
dominate the cosmic ray budget at all angular distances, having many UHECRs clustered 
near them. We discuss the likelihood of production of UHECRs in each of these sources. 

As pointed out before by different authors, the nearby radio galaxy Ce ntaurus A is an ex- 
cellent candidate for the acceleration of UHECRs in its jets and radio lo 



2009 



Hardcastle et all l2009l : iHondal l2009l : iRieeer fc AharonianI 2009 



jes (Moskalenko et al. 



events flGorbunov et al.ll2008l : iNagar fc Matuhchll2008l : Moskalenko et al.ll2009f ). We can test 



Dermer &: Razzaaue 



2OIOI ) . Cen A is a prominent 7-r ay source ( Abdo et al.ll2010al Jbl: iDermer fc Razzaque 



2010 : lPiranlL„_„. ^ , 

20101 ) , being associated with two UHECRs (lAbraham et al.ll2008l) and potentia llv even more 



the impact of Centaurus A in our results by removing it from the ILAC data and recalcu- 
lating the correlation signal of nearby ILAC AGNs without Cen A. The dashed line in the 
lower panel of Fig. [7] corresponds to the correlation signal without Cen A, which differs from 
the signal including Cen A (solid line in Fig. [7]) in the following way: the strength of the 
correlation is considerably weakened on separations smaller than ^ 5°, the minimum of the 
probability signal at 1^2 essentially does not change its associated probability and the min- 
imum at ip3 has its significance slightly decreased to 5cr. Therefore, statistically speaking, 
the impact of Cen A on the correlation is mainly in the significance of the short-distance 
associations. 

NGC 4945 is classified as either a starburst galaxy or an obscured Seyfert 2 nucleus 



showing no signs of the presence of relativistic outfiows (e.g., ISpoon et al.ll2000l : IChou et al. 
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20071 ). while ESO 323-G77 is a Seyfert 1 galaxy tha t displays evidence for the presence of 
a mildly relativistic outflow in its X-ray spectrum ( iJimenez-Bailon et al.ll2008l ). We have 
to be careful when associating these galaxies with the nearby PAO e vents, since neither of 



them display powerful outflows o r are energetically impressive (but see lBoldt &: Ghoshlll999 



Pe'Er. Murase fc Meszarod 120091 ) . Both galaxies are located quite close to Cen A and any 
events associated with Cen A could also be mistakenly attributed to NGC 4945 and ESO 
323-G77 in a correlation analysis, given their proximity (i.e. a case of accidental correlation). 

A better understanding of the nature of the 7-ray emission in NGC 4945 and ESO 323- 
G77 could shed light on the issue of their likelihood as accelerators of energetic cosmic rays. 
Given the available information on these AGNs, our analysis suggests that Centaurus A is 
the primary candidate for the production site of several UHECRs observed in the direction 
of its group of galaxies. We note that ESO 323-G77 shows evidence of the presence of a 
relativistically broadened Fe Ka ( jJimenez-Bailon et al.ll2008l ). which suggests that the su- 
permassive black hole in this gala xy could be rapidly spinning, with the Kerr black hole being 
a pot ential UHECR accelerator (JBoldt fc Ghoshlll999t Ivan PuttenI l2008t iDermer fc Menon 
2009h . 



5.2. Angular separations vs. deflection angles 

The correlation signal obtained in this work will be helpful to constrain the results 
of simul ations of the propagation of energetic cosmic rays origiri ating from extragalactic 
sources (JAbraham et al.l l2008t iLimal 120091 : iRyu. Das fc Kangl I2OIOI ) . One should be careful 
when mapping the separation angle ip that minimizes P for a set of reference sources into 
the actual deflection angles 9 suffered by the cosmic rays on their way towards Earth with 
respect to their actual sources. For instance, if the ILAC GZK sources that we studied are 
tracing the production sites of the observed PAO events, which is a tempting possibility 
given the high significances of association that we obtain, the double-peaked significance 
pattern in the correlation between ILAC GZK sources and PAO events raises a challenging 
question: Which of the angular scales V'l ~ 2°, '?/^2 ~ 7° or ■j/'s ~ 17° actually corresponds to 
a deflection angle? The answer to this question is beyond the scope of this work, but it is 
worthwhile to compare our results with numerical simulations. 



Ryu. Das fc Kangl (I2OIOI ) modeled the propagation of UHECRs from source AGNs to 



observers through the intergalactic magnetic field, using a model universe based on results 
from cosmological structure formation simulations. Interestingly enough, they obtained the 
mean value of deflection angle between the arrival directions of UHECRs and the actual 
sources of {9) ~ 14° — 17.5°. Therefore, our inferred angular distance ip^ ~ 17° - which cor- 
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responds to the maximum significance S{tp3) = bAa (see Figure [7 1) - is w ithin the theoretical 
range of mean deflection angles calculated by lRyu. Das fc Kang (120101 )! 



On the other hand, the separation angle between the arriv al directions of UHECRs a nd 
the position of the nearest "mock" AGN in the simulation of iRyu. Das fc Kangl (120101 ) is 
(V'sim) ~ 3.5° — 4°, therefore {9) ^ {ipsim)- Quite interestingly, when comparing the position 
of PAO events against those of all ILAC AGNs (left panel of Fig. ^ we find a much smaller 
separation angle that minimizes the probability of the observed configuration being due to an 
isotropic UHECR flux, ip ?a 2. 4°. Taken together , our r esults seem to be in rough agreement 
with the theoretical results of iRyu. Das &: Kang (120101 ). 



5.3. Relation to the supergalactic plane 

The SGP traces the distribution of matter in the local universe, with sey eral nearby 



Lahav et al.ll2000r) . as well as 



20081). For instance, 13 out of 



galaxy clusters and AGNs being concent rated near this plane 

the arrival directions of many UHECRs ( Abraham et al.ll2007l . 

the 27 (48%) PAO events are located within a distance of 15° to the SGP, with the average 

distance of UHECRs to the SGP of ^ 24.5°. 

Of the 9 ILAC AGNs with z < 0.048 in the PAO field of view, the 7 objects that 
correlate with energetic PAO cosmic rays within ip^ = 16.9° are located at an average 
distance to the SGP of (ip^gp) = 16.8°. For comparison, PAO events associated with the 
ILAC GZK AGNs are located at the average distance {'ipf^^) = 19.7°. This suggests that 
the nearby 7-ray bright AGNs associated with UHECRs on average fall very close to the 
nearby local overdensity of matter in the universe, even though the general distribution of 
ILAC AGNs does not follow noticeably the SGP. 



6. Concluding remarks 

We analysed the correlation between the positions of the 7-ray sources in the Fermi LAT 
First Source Catalog (IFGL) and the First LAT AGN Catalog (ILAC) produced during the 
first eleven months of operation of the Fermi Gamma-ray Space Telescope, and the arrival 
directions of the highest energy cosmic rays measured during 1.2 years of operation of the 
Pierre Auger Observatory. Using our cross-correlation test, we are particularly interested in 
the value of separation angle ip between the positions of PAO events and reference sources 
that maximizes the significance level S of rejecting the null hypothesis (that the observed 
configuration is due to chance assuming an isotropic flux of UHECRs). Our main results 
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from the cross-correlation analysis of each IFGL subset (in the PAO field of view) with the 
PAO events are summarized below: 

• When considering all 7-ray sources in the IFGL, ^max = 2.7cr at ip ^ 5.1° with less 
observed correlations than expected assuming the null hypothesis. 

• Galactic IFGL sources are not significantly correlated with the arrival directions of 
energetic cosmic rays {S < la). 

• Considering the ILAC AGNs we find S^ax = 2. 60- at ^ ^ 2.4°. 

• We find no significant correlation with ILAC blazars {S < la). 

• Considering the subset consisting of unidentified AGNs and non-blazars in the ILAC, 
5max = 3.5CT at ^ ^ 3.2°. 

• We find a significant association with ILAC sources closer than 200 Mpc, with the 
significance peaking at the separations ipi = 2.3° (4cr), iIj2 = 6.5° (5cr) and ip^ = 16.9° 
(5.4a). 

For the subset of nine ILAC AGNs closer than 200 Mpc {z < 0.048, roughly the distance 
to the GZK horizon), we reject with the highest significance the null hypothesis compared 
to the other subsets of IFGL sources (P^in = 5 x 10"^ at ^3 = 16.9°). Of the 27 PAO 
events, 17 are located within 16.9° of 7 nearby ILAC AGNs: Centaurus A, NGC 4945, ESO 
323-G77, 4C+04.77, NGC 1218, RX J0008.0+1450 and NGC 253. These "7-ray loud" AGNs 
are likely candidates for the production sites of UHECRs. 

The considerable separation angle that maximizes S could suggest that the trajectories 
of the particles on their way to Earth are significan tly defiected by the intervening magnetic 



fields, consistent with the numerical simulations of iRyu. Das fc Kand (120101 ). This may set 
important constraints for the understanding of the structure of magnetic fields along the line 
of sight. 

Given the small number of currently detected UHECR events, this work represents a first 
step towards understanding the relation among some of the highest-energy phenomena in the 
universe detected by PAO and Fermi. The data set that the PAO and Fermi will gather in 
the next few years should improve the correlation statistics of our work. Furthermore, cosmic 
ray observatories in the northern hemisphere such as the anticipated northern PAO site will 
offer a way of further identifying the sources of energetic cosmic rays. Thanks to PAO and 
Fermi, new perspectives to understand the relation among the highest-energy phenomena in 
the universe are being made possible! 
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Fig. 1. — Aitoff projection of the celestial sphere in Galactic coordinates, showing the 27 
highest energy cosmic rays detected by the Pierre Auger Observatory and the 1451 IFGL 
7-ray sources observed with Fermi (1037 in the field of view of PAO). The dashed line 
represents the border of the of the field of view of PAO, corresponding to a zenith angle of 
60°. The dot-dashed line corresponds to the supergalactic plane. [See the electronic edition 
of the Journal for a color version of this figure.] 
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Fig. 2. — Correlation signal between the arrival directions of UHECRs and the AGNs with 
z < 0.017 within the field of view of PAO in the 12th edition of the Veron-Cetty and Veron 
catalog. Upper panel: cumulative number of cosmic ray events as a function of angular 
distance ip (CCF; black line). The average chance expectation with the 68% confidence 
intervals is represented by the shaded region. Lower panel: the solid line correspond to the 
probability P that the number of observed correlated events is due to chance correlations. 
The dashed lines indicate the probabilities corresponding to the different significance levels 
for the null hypothesis. 
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Fig. 3. — Correlation signal between all IFGL 7-ray sources and energetic cosmic ray events. 
Upper panel: cumulative number of cosmic ray events as a function of angular distance ip 
(CCF; black line). The average chance expectation with the 68% confidence intervals is 
represented by the shaded region. Lower panel: the solid line correspond to the probability 
P that the number of observed correlated events is due to chance correlations. The dashed 
lines indicate the probabilities corresponding to the different significance levels for the null 
hypothesis. 
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Fig. 4. — Same as Figure [3] for the IFGL Galactic sources. 
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Table 1. ILAC AGNs with z < 0.048 and the associated PAO events 



IFGL 


I 


b 


Class^ 


Name'' 


x 


ID° 


Af(6.5°)'' 


Af(16.9°) 


name 


















J1305.4-4928 


-54.7° 


13.3° 


Sy 2 


NGC 4945 


0.002 


1 


4 


7 


J1325. 6-4300 


-50.4° 


19.4° 


RG 


Gen A 


0.002 


2 


3 


7 


J1307.0-4030 


-54° 


22.4° 


Syl 


ESO 323-G77 


0.015 


3 


3 


6 


J0047.3-2512 


97.4° 


-88° 


SB 


NGC 253 


0.001 


4 





1 


J1517.8-2423 


-19.3° 


27.6° 


BLL 


Ap Lib 


0.048 


5 








J0308. 3+0403 


174.9° 


-44.5° 


BLL 


NGC 1218 


0.029 


6 





2 


J2204.6+0442 


64.7° 


-38.6° 


Sy 1/BLL 


4C+04.77 


0.027 


7 


2 


3 


J1230. 8+1223 


-76.2° 


74.5° 


RG 


M87 


0.004 


8 








J0008. 3+1452 


107.6° 


-46.7° 


Syl 


RX J0008.0+1450 


0.045 


9 





1 



^Meaning of acronyms: RG 
starburst galaxy. 



radio galaxy, BLL = BL Lacertae, Sy 1 = Seyfert 1, Sy 2 = Seyfert 2, SB 



''Name of identified or likely associated source. 

° Identification of each source numbered in Figure |5] 

"^Number of PAO events potentially associated with each Fermi source at the given separati( 




Fig. 5. — Aitoff projection of the celestial sphere in Galactic coordinates, showing the PAO 
events (blue circles) and the ILAC AGNs. The square symbols represent the AGNs with 
z < 0.048 (i.e. within the GZK horizon) in the PAO field of view with circles of radius 16.9° 
centered at each square, while the x symbols denote all other AGNs. The numbers can be 
used to identify the ILAC GZK AGNs in Table [H The dashed and dot-dashed lines follow 
the notation of Figure [H [See the electronic edition of the Journal for a color version of this 
figure.] 
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Fig. 6.— Same as Figure [3] for the 453 ILAC AGNs. 
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Fig. 7. — Same as Figure [3] for the 9 ILAC AGNs within the GZK horizon. In the lower 
panel, the solid and dashed lines correspond to the correlation signal including and not 
including Cen A, respectively. 



